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1 Introduction 



The integrability theory for 1 + 1 - dimensional homogeneous systems of hydrodynamic type of 
the form 

n 

U i,t = ^2 a ij{ v ) u j,x, i = l,...,n, (1.1) 

j=i 

where u = (m, . . . , u n ) was developed in [H [2]. In the case of 1 + 2 - dimensional homogeneous 
systems a definition of integrability based on the existence of hydrodynamic reductions was 
proposed in [3]. For examples and a classification of integrable cases see [Jj and references 
therein. 

In contrast to the homogeneous case, there is no satisfactory criteria of integrability for 
non-homogeneous hydrodynamic type systems of the form 

n 

Uj,t = Oij(u) Uj,x + bj(u), i = l,...,n. (1.2) 

On the other hand, the following example 

u t = vu x H — , v t = uv x H — (1.3) 

v — u u — V 

found in [6J indicates the existence of integrable systems (11.2j) having properties unusual for 
1 + 1-dimensional integrable models. In particular, system ( II. 3p has only a few local infinitesimal 
symmetries and therefore the symmetry approach to classification of integrable 1+1-dimensional 
systems [7] is not applicable for systems similar to (11. 3p . It turns out that this system possesses 
infinite non-commutative hierarchy of nonlocal symmetries explicitly depending on x and t. 
These symmetries look like symmetries of the so called Gibbons-Tsarev systems found in [5]. 

Notice that system ( II. 3p was also derived in another context, in [Hi |9] where conditions for 
two quadratic Hamiltonians to be in involution were studied. 

It was shown in [10] that the general solution of (II. 3p can be described in terms of a 
conformal mapping of a slit domain to the upper half plane (cf. [TT1 IT2]). A class of special 
solutions was constructed in [T3]. All these results are related to the existence of the following 
Lax representation 

A — u — v 1 . 

(A — u)(X — v) {X — u)(X — v) 

for (II .3p . Here A is a spectral parameter. 

In this paper we consider systems (11.21) having Lax representations of the form 

tt s = /(u,A)#A, * t = ^(u,A)* A . (1.5) 
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For brevity, we call such systems integrable. Lax pairs similar to (11. 5p were considered in [13] . 
Note that Lax representations with A-derivatives of the \I/-function are typical for Painleve type 
equations. Possibly the technique developed in [16] for solving of dispersionless systems can be 
adopted to systems with such Lax pairs. 

To generalize example (II. 3p . (II. 4p we assume (cf. [T5]) that both functions / and g in (II. 5p 
have simple poles at A = Ux,...,u n . In general, we do not restrict ourselves to functions /, g 
rational in A but all examples in this paper are rational. 

In Section 2 we show that for any system (11.21) with a Lax representation of the kind 
described above the matrix a = (ay(u)) is diagonal and weakly non-linear [17]. The latter 
property means that the functions aa(u) do not depend on Uj. We establish a partial separa- 
tion of variables for functions / and g in ( 11. 5ft and write down a general functional equation 
describing these functions. 

In Section 3 we construct a large class of examples of integrable systems of the form ( 11. 2ft 
with arbitrary n. These systems depend on parameters. The main continuous parameters are 
the roots A, of a polynomial with constant coefficients. To each Aj a non-negative integer fcj 
is attached. In Section 3.1 the case is considered when all the k{ are zero. In Section 3.2 we 
describe a limit of formulas from Section 3.1 when some of the Aj coincide. In Section 3.3 we 
generalize results obtained in Section 3.1 to the case of arbitrary fcj. 

Section 4 is devoted to the case n — 2. We find all possible pairs of functions / and g 
appearing in the Lax representation ( 11. 5ft . Even for n = 2 the functional equation for functions 
/ and g is highly non-trivial. It turns out that it is almost equivalent to a functional equation 
describing the functions F, H in the fieldless Gibbons-Tsarev type systems 

diij = F(&, £,■) diU, didju = Q diudju, i ^ j, i,j = 1, N. (1.6) 

Here u, are functions of variables r 1 , . . . ,r N , and di = Such systems play a crucial role 
[3]II] in the integrability theory for 1 + 2-dimensional homogeneous hydrodynamic type systems. 
As an additional result we find new examples of Gibbons-Tsarev type systems (ll.6p . From our 
classification it follows that in the case n = 2 all systems with Lax representations (I1.5P can be 
constructed by the approach developed in Section 3. 

In Section 5 we present some explicit examples of integrable systems of the form (II .2p in 
the case n = 3. 

2 Lax pairs with movable singularities 

Our goal is a generalization of example (ll.3p . ( 11. 4p . In this example both functions / and g in 
(11.51) have simple poles at A = u and X = v. In the general case let us assume (cf. pj)]) that 
the functions / and g in ( 11.51) have simple poles at A = u\, u n . 

Computing {^ x )t and {^t)x by virtue of ( II. 2p and (II. 5p and equating coefficients at Uj )X , j = 
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1, ...,n, we obtain the following compatibility conditions: 

— = ^a Jl — , , = !,...,„, (2.7) 

<S^,df dg df 

Computing the singular part of ( 12. 7ft at A = itj, where j 7^ i, we obtain a^j = 5jjaj. We assume 
that the functions di are pairwise distinct. Now our system (jl.2p reads: 

Ui it = ai(u)u i)X + bi(vL), i = l,...,n. (2.9) 

The compatibility condition (12. 7p becomes 

R~ =(Xi TT' ^ = 1 > •••,«<• (2-10) 

Let 

f = T ^— + 0(\-u i ), g = -^— + o(X-Ui), 0,^0, ^0 
for A — >• u%. Computing the singular parts of ( 12. 10[) . we obtain 

4>i = diipi, 



dui dui 

for all i, j = 1, ...,n. Substituting the first of these equations into the second one, we obtain 
{djipj) Ui = cn4>j,Ui- From here it follows that if % = j, then 

da* 

— - = 0, 

dm 



otherwise 



dtp* 1 ddj , . 



dui di — dj du 



The first of these equations shows that the homogeneous system 

Uij = diU i>a; , i = l,...,n (2.11) 

is weakly non-linear. The compatibility conditions for the second equation means that this 
homogeneous system is semi-Hamiltonian. Weakly non-linear semi-Hamiltonian systems (12. lip 
were studied by E. Ferapontov in jTTJ . In particular, he classified these systems finding possible 
coefficients di in a closed form in terms of arbitrary functions of one variable. Namely, for any 
weakly non-linear semi-Hamiltonian system (I2.1ip the coefficient a, has the form: 



Here and in the sequel 



A 



/ 1 

<?i,iOi) 



\ g n -i,i(wi) 



where arbitrary functions of one variable 



1 \ 

<?l,n(«n) 



(2.13) 



• q n -l,n(u n ) ) 

The meaning of the indexes on A in (2.12) is 



as follows. For any matrix M we denote by My its minor. In other words, My is obtained 
from M by deleting its i— th row and j— th column. By Mj we denote the matrix obtained from 
M by deleting its z— th row. 

For functions (12. 12ft the general solution of (I2.10p is given by 



./ 



detP 
det A : 



where 



P 



hi(X,ui) ... h n (X,u T 
A„-i 



Q 



detQ 
det A' 

/ii(A,«i) 



(2.14) 



h n (X,u r 



Here hi, h n are arbitrary functions of two variables. Notice that P n = Q n . According to our 
assumption each function hi(X,u) has a simple pole at A = u. Substituting these expressions 
for /, g into (12. 8p and computing the singular part at A = Mj, we get 



det Qn,i\\=Ui 

det A n _i i 



i — 1, n. 



(2.15) 



Substituting these into (I2.8p . we obtain a functional equation, which can be written in the 
following form: 



det 



hi(X,Ui) x ... h n (X,u n ) x 

Qn 



J +^(-l)^ 1 ^(A,u i ) Ui detQ n , J | A=Ui + (2.16) 



1 



det A 



E (-l) fcl+fc2 detQ n , fc2 | A=Ufc2 det( 



/ii(A,tii) 



h n (X,u r 



^fci+i 



l<fcl<n-l,l<fc 2 <n 

Remark 1. Expanding numerators of (12.141) by the first row, we get 

/ = /ii(A,«i)0i + ... + h n (X,u n )(p n , g = /ii(A,Mi)^i + ... + h n (X,u n )ij. 

where 

\ . , ■ 



-jVi-l ^n,i 

A n 



0i = (- 1 ) "a — S 

Note that according to the Cramer rule, 0j, ^ are the solutions of the following system of 
linear equations: 



q j>X (j)X + ... + q jtn <f) n = 8 n -2,j, Qj,l1pl + ••• + qj,ni>n = 8 n -l,j, 



(2.17) 
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where j = 0, ...,n — 1, %,i{u) = 1 and 8ij is the Kronecker delta. □ 

Remark 2. For all known examples the functions hi(X,u), qj,i{u) don't depend on i. We 
conjecture that this property holds for any solution of (12.161) . This conjecture is proved in the 
cases n = 2, 3. □ 



3 A class of integrable systems with arbitrary number 
of components 

The compatibility condition for (jl.5p read^] 



dt dX dx d\ 



(3.18) 



In this section we assume that / and g are rational functions in A. In addition, we suppose 
that both / and g are divisible by a polynomial S(X) with constant coefficients. At first glance 
this is not true for our basic example (II. 4p . Nevertheless, under the transformation 



A ->• A 



aX + b 



(3.19) 



cA + d 

with constant a, b, c, d a common multiplier (cA + d) 3 in / and g appears. 

One can construct conservation laws for a system of the form (11.21) using such a Lax repre- 
sentation (ll.5p as follows: 

Lemma 1. Suppose that / = S(X)f\, g = S(X)gi, where S(X) is a polynomial with 



constant coefficients and with roots Aj, X p . Set \ii 

dx 



f: 



1 A=A, 



andui = g 1 \\=\ i . Then 



dt 



Moreover, if A, is a root of multiplicity k and /ijj 
0, k — 1, then 



V, 



&9x 
dXi 



A=Ai 



d/i 



dv, 



1,3 



dt dx 

Proof. Substituting our expressions for /, g into (13.18j) . we get 



df 



^ + 5(A)^A = ^ + 5(A)— 



dgi 



(3.20) 
where j = 
(3.21) 

(3.22) 



To complete the proof we set A = Aj here. In the case of multiple roots we compute derivatives 
of (I3.22p with respect to A and then set A = Aj. □ 



1 Note that this equation can be written as |£ — ^ + [/, g] = 0, where [f,g] = /Jf — g|f is the bracket in 
an algebra of vector fields. Therefore p. 181) can be considered as zero curvature equation for a connection in an 
Di f /^-bundle, where Diff 1 is the group of diffcomorphisms of one-dimensional manifold with coordinate A. 
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3.1 Principal series of examples 



Let 

_ S(\)P(X) _ S(\)Q(X) 

f ~ R(X) ' 9 ~ R{X) ' 

where S, P, Q, R are polynomials in A and all coefficients of S are constant. Then equation 
(I3.18P reads as 

^f + ^-^ + W^-^O. (3,4, 



dt dt dx dx \ d\ d\ 

Proposition 1. Let 

S(X) 

hi(X,u) = - _ , qj,i{u) = u 3 

if i = 1 , .. . , n, j = 1 , . . ., n — m — 3 and 

n+l— m 

Cfc 



g„_ m _ 3+iii (u) = ^ \ (3.25) 
fc=i 

if i = l,...,n, j = l,...,m + 2. Here m is a fixed integer such that < m < S(X) = 
(A — Ai)...(A — A n+ i_ m ), where Aj are pairwise distinct constants, and cy are arbitrary constants 
such that the matrix (cy) has rank m + 2. Then formulas (12.131) - (12.151) define a system (12. 9p 
possessing a Lax representation (11.51) . 

Proof. Expanding the determinants in the numerators of (12.141) by the first row, we get 

/ = W^ + ... + ^U, 9 = S(A) (*_ + ...+ *!_), (3.26) 

where (pi, 4> n ,ipi, ...,ip n satisfy the linear system (I2.17P (see Remark 1 in Section 2). The 
part of this system corresponding toj = 0, n — m — 3 has the form 

+ ... + <0 n = 0, tt}Vi + - + <i/> n = 0. (3.27) 

From these equations it follows that 

S(X)P(X) S(X)Q(X) 
(X-ui)...(X-u n ) (A - ui)...(X - u n ) 

where P, Q are polynomials in A of degree m+1. Indeed, using expressions (I3.26P and expanding 
^ and ^ with respect to powers of A -1 , we get left hand sides of (I3.27P as coefficients at A - - 7 , 
j = 1, ...,n — m — 2. The left hand side of (I3.24p is a polynomial in A of degree n + m + 1. 
Equating its coefficients to zero, we obtain n + m + 2 non-homogeneous differential equations 
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of hydro dynamic type for u\, ...,u n . To obtain a system of the form (jl.2p we have to prove that 
there exist m + 2 linear dependence relations between these equations. Let 



fJLi 



Xi - Mi 



+ ... + 



A,- - u. 



V>1 , , 



Aj - Mi 



Aj - u„ ' 



(3.29) 



where % = 1, n + 1 — m. According to Lemma 1, we have n + 1 — m differential equations 
(I3.20p . On the other hand, the relations 



C 1J U 1 + ••• + C n _|_i_ m .j/! n +i- 



c l,j v l + ■•• + Cn+l-mj^n+l-r 



m+2,j 



(3.30) 



for j = 1, ...,m + 2, are fulfilled. Indeed, substituting the expressions ( 13. 29ft for u i; z/, into 
(I3.30p . we get the part of the linear system (12.171) corresponding to j = n — m — 2, n — 1, 
where the functions q^i are given by (I3.25p . Relations (I3.30p give us m + 2 linear dependence 
relations between equations ( 13. 20 p . namely 

Cljiil^^t — (^lja;) + ••• + Cn+l-m,j(\l J 'n+l-m)t ~ {^n+l-m)x) = 0. 



Since (fii)t — {vi) x is proportional to the left hand side of ( 13.241) . where A is set to Aj, we obtain 
m + 2 linear dependences between coefficients of the left hand side of (13.241) . Note that Uj (resp. 
Vj) are given by the formula (I2.14p for / (resp. for g), where hj(X,Uj) are replaced by x .\. D- 

Here we present a more explicit form of the functions / and g. Let S(X) = (A — Ai)...(A — 
A n +i-m)j where Aj are pairwise distinct constants and R(X) = (A — «i)...(A — u n ). Functions 
/, g can be written as 

f(ui, u n+ i_ m , t>i, v m , A) = (3.31) 



5(A) Ei<n 



<...<j m +i<n+l— m V*H— 



i(A - X h )...(X - X im+1 )R(X h ) 1 ...R(X im+1 ) 



-i 



R W ^2l<h<...<i m+2 <n+l-m A ii...W2- R (Aii)" 1 ---R(A im+2 )" 1 

g(ui, u n+ i_ m , v m , A) = 
'S'(A) Y J i<i 1< ... < i m+1 <n+i-m^h-i m +i(^ ~ X h )...(X - A im+1 ) J R(A il )" 1 ... J R(A Wl )~ 1 



<...<i rri _l_2 



_i_ m A il ... im+2 _R(A il ) 1 ...R(X im+2 ) 1 



where 



1.4+1 = FI ( A ^-A i(3 )det 

l<a</3<m+l 



^ii.-W = II ( A ^-A i/3 )det 

l<a</3<m+l 



Ci u l 



tt «m+l \ 
C «m + l,l 



c i m+ l,m / 

c Wi>i 

•• c i m+ i,m / 
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A, 



<,...i m+2 = I I (K a - X i0 )det 

l<Q</3<m+2 



/ 


Oil 


a im + 2 


\ 














C »m+2,1 




\ 


Qi,m 


• C »m+2,m 


/ 



and a,i, bi, Aj, q,- are constants. In particular, if m = 0, we have 



/(ui, ...,u n+h A) = 
gr(u l3 ...,u n+ i, A) 



5(A)Er=i 1 (A-A 4 )6, J R(A t )- 1 



(3.32) 



Note that this form is useful for changing coordinates. For example, we can choose R(\) 
v 1 +v 2 X + ... + v n X n ~ l + X n . 



3.2 Degenerations 

The construction of Proposition 1 admits a limit when some of Aj coincide. All formulas are 
valid in this case except (I3.25p . The formula (I3.25P should be rewritten as 



,n—m 



Qn—m—3+j,i (U) = \ v , 7 \ , (3-33) 

(U - Al)...(lt - Xn+l-m) 

where cy are constants. In (I3.33j) some of Aj may coincide. An analog of (I3.25|) can be 
obtained from ( 13.331) by the partial fraction expansion of rational functions in u. Note that if 
Xi = Aj+i = ... = Xi + k, then the corresponding densities of conservation laws (13.2 ip read 

= T\ + •- + 7a 17^7' v hi = T\ 7T77 + - + 



where j = 1, k + 1. 

It is clear from ( II. 5p that functions /, g of the form (I3.23P with degrees of polynomials 
P, Q, R, S equal m + 1, m + 1, n, n + m respectively transform as 

, , (cA + rf) 2 5(A)P(A) (cA + d) 2 S(A)Q(A) 
/_ > / = fl(A) ' ^ = W) 

under transformations (13 . 1 9[) . In particular, one of the values of Aj can be sent to infinity. If 
Aj = oo and this root has multiplicity k, say A n+ i_ m = ... = X n+ 2- m -k = oo, then 5(A) = 
(A — Ai)...(A — A n+ i_ m _fc) and (I3.33P takes the form 

— I — l l — 71 — 771 

( \ _ Cl 'J ~"~ C 2,i U i" ••• ' C n+l-m,jU 



u - Xi)...(u - A n+ i_ m _ fc ) 
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Consider the case when all Ai = ... = A rt+ i_ m = oo. We have S(\) = 1 and 

Qn-m-3+j,i( u ) = + ^2,jU + ... + C n+ i_ mJ « n "\ 

Since we can replace the equations in system (12.171) by any of their linear combinations, in 
particular we can subtract linear combinations of equations (I3.27j) from other equations, we 
can assume without loss of generality that 

q„_ 3+j>i (u) = CiX~ m ~ 2 + c 2j u n ~ m - 1 + c 3j u n ~ m . 

Therefore we can only have m = or m = 1. If m = and our constants Cy are chosen in such 
a way that 

q n -2,i{u) = u n ~ 2 , = u n ~ l , 

then our system takes the form 

Ui,t = 22 u j u i,x + Y\_( Ui - * = 1; 

This system has appeared in [Sj, where a different problem was studied. If n = 2, then this 
system coincides with (II. 31) . 

3.3 General scheme 

The main idea of our construction from Proposition 1 can be described as follows. Let L(\) be 
the difference between the left and right hand sides of ( 13.181) . If /, g are given by (13.231) . then 
L is a rational function in A and its numerator is the left hand side of ( 13.241) . If the degrees of 
the polynomials P, Q, R, S are m + 1, m + 1, n, n + 1 — m respectively, then the degree of 
the numerator of L equals n + m + 1 and therefore the identity L = is equivalent to a system 
of n + m + 2 non- homogeneous hydrodynamic type equations. To provide m + 2 linear relations 
between them we impose constraints of the form 

cijL(Ai) + ... + c„ + i_ TO jL(A„ + i_ m ) =0, j = 1, m + 2, 

where the cy are constants. Since L(Aj) = f(\)t ~ g(\)x, these constraints follow from 

c i,j/(Ai) + ••• + c n+ i_ mj /(A n+ i_ m ) = dj, Cijg(\i) + ... + c n+ i_ m j^(A n+ i_ m ) = bj, (3.34) 

where j = l,...,m + 2 and dj,bj are arbitrary constants!! The coefficients of P and Q are 
uniquely expressed from (13.341) in term of coefficients of R. 

It is possible to generalize the construction described in Section 3.1 in the following way. 
We still assume that / and g are given by (I3.23p . where the degrees of polynomials P, Q, R, S 

2 Without loss of generality we can set ctj = S m +i,j, bj = § m ,j- 
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are m + 1 , m + 1 , n, n+l — m and < m < n — 1 . In the case m<n-lwe use representation 
(13T26|) . If m = n - 1, then 

/ = S(A) (fo + + ... + j^-) , , = S(A)(*+ *+...+ *!_Y (3.35) 

V A — Ui X-u n J \ A- Mi X-u n J 

In our generalization we suppose that (13. 34ft is valid for j = 1, ...,m + 2 — k. We look for the 
remaining k linear constraints in the form 

di t j,if(\i) + d2,j,if'(\i) + ••• + di t j t i ~ (/_ 1 )t'' > = 0, 

i-i ( 3 - 36 ) 
di,j,ig(K) + d2 t j l ig'(\i) + ... + (;_i)t' = 0) 

where j = 1, k iy i — 1, n+ 1 — m and fcj + ... + fc„+i_ m = fc. The coefficients d s _j^ have to 
be found from the identity 

dt^LiXi) + dwL'iXi) + ... + diJ,i L l^f = 0- ( 3 - 37 ) 

Let 

/ = (A - Xi)f + (A - A,) 2 /i + 9=(X- \)9o + (A - Xifgi + 

L = L + (A - A»)Li + (A - A,) 2 L 2 + ... . 

We omit the index i in the coefficients of these Taylor expansions for simplicity. Substituting 
the Taylor expansions into (13 . 1 8j) . we obtain 

Lq = fo,t — go,x-, 

Li = fi, t - g%, x + ^2 (* _ 2k )(fkgi-k - fi-k9k), i > o. 

Substituting these expressions into (I3.37p . we see that the terms fj )t , gj ;X cancel out by virtue 
of (I3.36p . The remaining parts of the expressions are bilinear in fj, gj. The vanishing of such 
parts leads to constraints for the coefficients d s ^ t i. Given a set of coefficients d SJj j satisfying 
these constraints, we substitute expressions (13.261) for /, g into (13.361) and obtain a set of 
k — ki + ... + k n+ i_ m linear equations for fa, We combine this system of linear equations 
with (I3.27P and (I3.30p . where are given by (13.291) and j = 1, ...,m + 2 — k. The whole 

set of linear relations guarantees that the identity L = is equivalent to a system of the form 

(O). 

We do not describe here all admissible sets of coefficients d s ,j,j in ( 13.36P but just discuss two 
admissible cases. 

Case 1. Suppose that (I3.36P takes a form f s = d Sji f , g s = d Sti g , s = 1, Then we 

have also L s = d s> iLo, s = 1, ki. 
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Let m < n — 1. In this case /, g are given by (I3.26p . where <j>i, <j> n , if)i, ipn are defined 
as the solution of the linear system combined by the following three parts: 

Part 1. The system (I3.27P with j = 0, n — m — 3. This system is empty if m = n — 2. 

Part 2. The system ( I3.30p . where //j, ^ are given by ( 13. 29ft and j = 1, ...,m + 2 — k. 

Part 3. The system 

' Hl ^ + - + ( 7~\ 1 ^ - ^ = 0, (3.38) 



1 p * v+-+ l 1 v^ -t^*-k=o. 



(A, - Xi-uxJ \{Xi~ u n y +1 Xi - u r 

where j = 1, k iy i — 1, n + 1 — m and fci + ... + k n+1 _ m = k. 

If m = n — 1 then we have to take = n — 1. The functions /, g are given by ( 13.351) . where 
(fio, (fi n ,ip , ...,ip n is the solution of the linear system combined by the following two parts: 

Part 1. The system (I3.30p . where m, z/j are given by 



^ = 0O + V~ \- - + T— ^ > ^ = ^0+, 

Aj - Mi Aj - W n Aj - Mi Aj - M n 

and j = 1, 2. 

Part 2. The system 

- Pj<i<h + ( 77 ^7TT " — 1 0i + - + ( 77 ^7TT " — 1 = 0> ( 3 - 39 ) 

V(^- m i) j Xi-uxj \(Xi-u n y +i Xi-u r j 

where j = 1, hi, i — 1, 2 and k\ + k2 = n — 1. □ 
Case 2. Suppose that (13.361) takes the form 

f s = a s ,ifo + b s ,ifii 9s = a a ,i9o + b a ,i9u s = 2, (3.40) 
where the coefficients ajj and 6 Pj9 satisfy the following equations 

(s - 2)a s _ M + 6 Sji + >^ (s - 2k)(b k ,ia s -k,i - b s - k ,iak,i) = 0, s = 2,3, .... 

2<fc<f 

In particular, 62.1 = 0. Then we have also L s = a Sj jLo + b Sj iLi, s = 2, fcj. The explicit form of 
a linear system for l5 <p n , ipi, ^> n , where m < n — 1 (resp. for O , <f) n , ip , where 
m = n — 1) can be obtained straightforwardly by substitution of ( I3.26P (resp. ( 13.35(1 ) into ( 13.4=0(1 
and combining these equations with (I3.27P and ( 13.30(1 . □ 
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Remark 3. For any admissible coefficients d s j,i, if Aj are distinct, then we must have 

2m + 1 — n < k < m, k\ + ... + /c„+i_ m = k. (3-41) 

Indeed, we must have at least two equations of the form (13.301) for 0; (resp. for so k < m. 
On the other hand, the rank of the matrix (cj^) should be equal to m + 2 — k. Therefore, 
m + 2 — k<n + l — moT 2m + 1 < n + k. In particular, if m — n — 1, then k = n — 1. □ 

Consider the case n = 2. It is easy to see that (I3.4ip admits only two solutions: 

1. m = 0, k = k\ = k 2 = &3 = 0. This is the case of Proposition 1. Explicit formulas are 
given by (I4.46p . 

2. m = 1, k — 1. Without loss of generality we can set fci = 1, &2 — 0- There are only two 
possibilities for coefficients d s j^ and both are described by the cases 1 and 2 above. In case 1 
we have fi = afo, g\ = ago where a is an arbitrary constant and fi, gi are coefficients of the 
Taylor expansion of /, g respectively at A = Ai- Explicit formulas are given by f )4.49p . In case 
2 we have f 2 = af , g 2 = ag , where a is an arbitrary constant and fi, g^ are coefficients of the 
Taylor expansion of /, g respectively at A = Ai. Explicit formulas are given by (I4.47p . 



4 Classification in the case n = 2 

Consider the case n = 2. In this case we do not assume that / and g are rational functions in 
A. Denote U\ — u, u 2 = v. Expanding the functional equation (" 12 . 161) a neighborhood of the 
diagonal v — u, we obtain qi,i(u) = qi t 2(u) and hi(X,u) = h 2 (X,u). Let 

qi,i(u) = qi, 2 (u) = — hi(X,u) = h 2 (\,u) = h ^ ,7 ^\ 
a[u) a[u) 

Then system (11.21) takes the form 

u t = a{v)u x + h(u, v), v t = a{u)v x + h(v,u) (4.42) 
and the corresponding Lax pair ( II. 5ft is defined by 

h(X,u) - h(X,v) a(v)h(X,u) - a(u)h(X,v) 

f{u, A) = — — — , g(u, v, A) = — — . (4.43) 

a[u) — a{v) a[u) — a{v) 

Here h is a function with a simple pole on the diagonal. Formula (I2.16P yields the following 
functional equation for this function: 

h(X, v)h(X, u)\ — h(X, u)h(X, v)\ + h(u, v)h(X, u) u — h(v, u)h(X, v) v — 
h(u,v)a'(u) — h(v,u)a'(v) 



a{u) — a(v) 



■(h(X,u)-h(X,v)) = 0. (4.44) 
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Remark 4. Suppose (14.441) holds and h has a simple pole on the diagonal. Then the 
following Gibbons- Tsarev type system with two fields u, v [6j [3j 0] is compatible: 

diPj = h(pj,pi)diU, diV = a(pi)diU, didju = a(p) - a(p) —diudjU. (4.45) 

Here i ^ j = 1, ...,N; p\, ...,px,u are functions in ri, ...,r/v, d{ = and AT is arbitrary. Note 
that there exist Gibbons-Tsarev type systems that have a slightly different structure. □ 

Proposition 2. For each Gibbons-Tsarev type system of the form (14.451) one can construct 
a non-homogeneous system (12. 9p with n = 2 possessing a Lax representation. 

Proof. Set iV = 2 in (I4.45P and consider pi, p 2 as functions of u, v. We have 

diPj = (Pj)udiU + (pj) v diV = (pj)udiU + (pj) v a(pi)diU = h(p j} pi)diU 

or ipj)u J r{Pj)vO'{pi) = h(pj,pi), where i j = 1, 2. Moreover, it is known [1] that each Gibbons- 
Tsarev type system admits so-called dispersionless Lax operator. It is a function L(\, r±, r^) 
defined (up to transformations of the form L — > q(L)) by the following system 

diL = h(\,pi)L x diU, i = l,...,N. 

Note that this system is compatible by virtue of (14.451) . Taking iV = 2, we may consider 
L as a function of u, v, X. As above we obtain L u + L v a(pi) = h(\,pi)L\, i = 1,2 or 
L u = fipi,P2, A)L A , L v = g(pi,p 2 , A)L A , where 

_ a(p 2 )h(X,p 1 ) - a(p 1 )h(X,p 2 ) = fo(A,pi) ~ H\P2) D 

a(p 2 ) - a(pi) a(pi) - a(p 2 ) ' 

Remark 5. In [8] several systems of the form (14.421) possessing a conservation law with 
density and flux that depend on u,v,u x ,v x were found. The existence of such a conservation 
law was proposed as an indication of complete integrability. Most of the systems from [8] do 
not satisfy our functional equation (14 .44 p . □ 

Let us present several solutions of the functional equation (I4.44p . 

Example 1. The functions 

S 3 (u) S 3 (v) , . P 2 (u) , , 

h{u > v) = w 2 (v)(u-vy a{u) = w^Ty (446) 

where S3, P 2 and W 2 are arbitrary polynomials of degree 3, 2 and 2 correspondingly, satisfy 
(EUD - For equation ([L3D we have S 3 (u) = 1, W 2 (u) = -1, P 2 (u) = -u. □ 

Example 2. Let 

Q(u) Q{v)P{u) (p(v) 2 P{u) + {u- v)Q(vj) s(u) 

Ku, v) = )— , L a {u) = -±±, (4.47) 

R(v)(u — v) R(u) 

14 



where P(u) = p\U + po, Q{u) = q\u + go are arbitrary polynomials and the coefficients of 
polynomials S{u) = s 2 u 2 + S\U + s and R(u) = r 2 u 2 + r\U + r satisfy the same linear equation 

2p <?o^2 - (po9i + Pi?o)»i + 2pigix = 0. (4.48) 

Then the functions h, a satisfy (14.441) . Notice that (I4.48P means that the double ratio (pq, ziz 2 ) 
equals —1, where p, q, zi, z 2 are roots of the polynomials P(u), Q{u) and x 2 u 2 + x\U + xq re- 
spectively. □ 

Example 3. The functions 

Q{u) Q{v)R{u) (r{v)R{u) + k{u - vj) c (v) 

Hu, v) = -±y L a{u) = W (4.49) 

1 {v)(u — V) 1 (u) 

where Q(u) = q\u + go S(u) = S\u + Sq, R(u) = r\U + r , T{u) = t\U + to ar e arbitrary 
polynomials, satisfy ( I4.44j) . □ 

It is easy to verify that the classes of solutions described in Examples 1-3 are invariant with 
respect to the transformations 

^ _^ k x u + k 2 ^ _^ kiv + k 2 ^ ^ 

k 3 u + ki ' k 3 v + k 4 

and 

x —> r\x + r 2 t, t — > r 3 x + r 4 t, (4-51) 

where fc« and are arbitrary constants. 

Remark 6. Using transformations (I4.50p . (14.511) . we can bring some of polynomials from 
Examples 1-3 to a canonical form. For instance, in the generic case of Example 1 one can bring 
polynomial S 3 to S 3 (u) = u(u — 1) . □ 

Theorem 1. Any solution of (14.441) that has a simple pole on the diagonal is given by 
P~46j) , g33) or (1449]) up to transformations (14331 . 

Proof. Let 

h(z, x) = ^ziM + ao M _|_ ai(x)(z - x) + .... (4.52) 
z — X 

Relations (I4.42p . (14.431) . (14 . 44 j) admit arbitrary transformations of the form 

u -)■ <f)(u), v -)> (4.53) 

Normalizing a(w) by u with the help of (I4.53P and expanding our functional equation at A = u, 
we obtain 

h(u,v) u = — ^ r(a-i(u)h(v,u) + {{u - v)a'_ x {u) - a-i{u))h{u,v)\. (4.54) 

2a-i{u){u — v) \ / 
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For the function 

, . h(u,v) ,. rrN 

A/a_i(Mja_i(v) 



equation (I4.54p takes the form 



q(v,u) - q(u,v) 

q{u,v) u = r . 4.56 

2(u — v) 



From ( I4.56P it follows that q(u,v) u = q(v,u) v . Differentiating (I4.56P by v and eliminating 
q(v,u) and q(v,u) v , we arrive at the following Euler-Darboux equation 

3 q u 1 Qv i A r<-!\ 

q uv = - (4.57) 

2u — v 2u — v 

for q(u, v). 

We are interested in solutions of (I4.57P of the form q(u, v) = —— + G(u, v), where G(u, v) is 
holomorphic on the diagonal u = v. It is easy to verify that for any such a solution the function 
G is of the form 



2k-l 



for some function s(x). The functions a_i(x),ao(x) from (I4.44p and s(x) are related as follows 

a = —y- + sa_i. (4.58) 

We will show that for any solution (I4.52p of (I4.44p the corresponding function s is rational 
(in contrast with a_i,ao). For any rational s(x) the corresponding function G can be easily 
reconstructed in closed form (cf. [9]). In particular, if s(x) = l/(x—k) then G(u, v) = T(u, v, k), 
where 

T(u, v, k) = 2 =. (4.59) 

{s/u — k) + vt) — k)yv — k 

The multiple pole six) = l/(x — k) 2 corresponds to — and so on. 

ok 

Substituting (I4.52p . (I4.58P into (I4.44p and expanding in a small neighborhood of u — v — z, 
we find that all coefficients a^, with i > are uniquely determined through a_i and s. These 
two functions satisfy an overdetermined system of ODEs. Eliminating a_i, we arrive at another 
system for s only. This system contains two ODEs of fifth and fourth orders. Differentiating 
the fourth order ODE and eliminating the fifth derivative by fifth order equation, we get several 
more fourth order ODEs. The final system is so complicated that usual algorithms of computer 
algebra such as differential Groebner basis technique certainly do not work. Fortunately, the 
system admits a group of point symmetries. 
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Although we fix a(u) = u, the equation f |4.44j) still admits symmetries. Namely, after any 
transformation of the form (14.511) the function a(u) becomes fractional-linear and we can bring 
it back to u by an appropriate transformation (14.501) . The existence of this GL(2)-action implies 
the symmetry group 

(kill + k 2 



s[u) 



2(k 3 u + k 4 ) (k 3 u + k 4 ) 2 

of the ODE system for s(u). Taking the simplest differential invariants 



s" + \2ssf + 16s 3 s'" + 24ss" + 144s 2 s' + 144s 4 



x = 77, , > V 



s / + 2s 2)3/2 ' » (s' + 2s 2 ) 2 

of this action for new dependent and independent variables, we can reduce the order of equations 
by 3. The equations of fourth order turn into very complicated first order equations for the 
function y = G(x). Considering these equations as polynomials in x, G and G' and eliminating 
G' by using the resultant several times, we find that the whole system is equivalent to the 
following single algebraic equation 

G 6 + 288 G 5 - 216 (5x 2 - 124) G 4 + 54 (13x 4 - 944x 2 + 16768) G 3 

-3888 (27x 4 + 40x 2 - 3632) G 2 + 5832 (63x 6 - 1124x 4 + 1856x 2 + 17920) G 
-729 (343x 8 - 9376x 6 + 69120x 4 - 118784a; 2 - 409600) = 0. 
The general solution of this third order ODE for s(u) is given by 

/ x 1 1 kr, 

s(u) = :— + ^ + t V-^, 4.60 

A(u — ki) 4(w - k 2 ) (u - ki) 2 

where fcj are arbitrary constants. Besides ( I4.60P there is a special solution of the form 

s{u) = 1 + h. (4.61) 
4{u — ki) 

The next step is to find the function a_i from the initial ODE system for s(u) and a_i(w). 
It turns out that if s is given by (I4.60p with k 3 = then this system is equivalent to a 
single non-linear ODE of fourth order for a_i. It can be linearized by the substitution a_i = 
(u — ki)(u — k 2 )P(u) 2 . Solving the linear ODE for P, we find 



a-i(u) — (u— ki)(u - k 2 ) [Pi{u)\/u - k x + P 2 {u)y 'u - h 
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where Pj are arbitrary first degree polynomials. The function g(u, v) is reconstructed from s as 
follows | ^ ^ 

g(u, v) = h -T(w, v, fci) + -T(m, v, k 2 ), 

u — v 4 4 
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where T(u,v,k) is given by (14.591) . The function h(u,v) is defined by (I4.55p . The radicals 
\fu — k\ and \Ju — k 2 can be removed by an appropriate transformation of the form (I4.53P and 
as result we get the solution described in Example 1. 

Example 2 corresponds to (I4.60P with k 3 ^ 0. In this case we have 

( \ 1 , l rr( / \ , 1 m/ 7^,7 ST(u,v,ki) 
g{u, v) = h -1 (u, v, ki) + —1 (u, v, k 2 ) + k 3 - 



u — v 4 ' 4 ' " dk\ 

and a_i = a(u — ki) 3 (u — k 2 ) 2 , where a is arbitrary constant. 

The degeneration k 2 = ki, k 3 ^ of (I4.60p gives rise to Example 3. In this case a\ = 
(u — kx) 3 Pf, where Pi is an arbitrary first degree polynomial. 

Solutions of (I4.44p corresponding to other degenerations can be written in one of the forms 
f)4.46p . (I4.47P or (I4.49P with special polynomials there. In particular, the degeneration k 2 = k±, 
k 3 = corresponds to (14.46)) . where the polynomials P 2 and Q 2 have a common root. □ 



5 Examples in the case n = 3 

In the case N = 3 integrable systems have the form 

u t = ai(v, w)u x + bi(u, v, w), 

v t = a 2 (u,w)v x + b 2 (u,v,w), (5.62) 
w t = a 3 (u, v)w x + b 3 (u, v, w), 

where 

B{u) - B{v) 

Oi(u,u) = a 2 [u,v) = a 3 (u,v) = —— — -, 

A(u) — A(v) 

X(v,u)-X(w,u) , X(w,v)-X(u,v) 

bi(u, v, w) = — — - — , b 2 (u, v, w) - 



A(v)-A(w) ' ' v ' ' ' A(w)-A(u) ' 

X(u,w) -X(v,w) 

b 3 {u,v,w) = — — . 

A(u) — A[v) 

The Lax representation of the form (II. 5p is given by 

X{u, X)(A(v) - A{w)) + X(v, X)(A(w) - A{u)) + X(w, X)(A(u) - A{v)) 



f(u,v,w,X) 
g(u,v,w,X) -- 



B(u){A{v) - A(w)) + B{v){A{w) - A(u)) + B(w)(A(u) - A{v)) 
X(u, X)(B(v) - B(w)) + X(v, X)(B(w) - B{u)) + X(w, X)(B(u) - B(v)) 



B(u)(A(v) - A(w)) + B(v)(A(w) - A(u)) + B(w)(A(u) - A(v)) 
In the above formulas the functions X(u, A), A(u), B(u) satisfy a complicated functional equa- 
tion, which follows from (I2.16p . 

Remark 7. Taking into account transformations (14.51)) . we see that the following group of 
afline transformations is admissible: 

A -)■ dA + c 2 B + c 3 , B -)■ c 4 A + c 5 B + c 6 , 
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5.1 A class of solutions 



Let us describe explicitly all integrable equations with the function X of the form 

A — u 

for some functions R, S. In this section we do not assume that the functions R and S are 
polynomials. Note that the transformation of X under (I4.50p is given by 



R(u) -+ (k 3 u + h)R(^^), S{\) -+ (k 3 X + hfs(^±^). (5.64) 



It can be straightforwardly verified that all examples found below belong to the class of equa- 
tions described in Section 3. 

From (12.161) it follows that the functions R, S satisfy several ODEs linear in S. The simplest 
of them are 

6RR"S® + 5{3R'R" - RR"')S {4) = 0, 
R 2 S {6) + URR'S® + 15(2R' 2 - RR")S { ^ = 0, (5.65) 
RR"S® - 8RR'"S® + 5(RR"" - 2R'R"' - 3R" 2 )S {4) = 0, 

and 

i? 4 + 5R 3 R' - 20R 3 R" + 10R 2 (3R'R" - RR'") S" - (5.66) 
5R(R 2 R {i) - 12RR" 2 + 12R' 2 R") S' - {R 3 R^ - 30R 2 R"R"' + 90RR'R" 2 - 60R' 3 R") S. 

There are two different possibilities: Case A: S"" = and Case B: S"" ^ 0. In Case 
B the determinant of the system of linear equations (I5.65P for should be zero. 

This leads to a fourth order ODE for R, whose solution is R = W2/W1, where W2 and W± are 
arbitrary polynomials of degree 2 and 1, correspondingly. Using a transformation of the form 
(15.64p . we bring W x to 1. 

In Case A besides (15361) and = we use one more equation linear in S. This equation of 
order 3 in S and order 4 in R is rather complicated and we do not present it here. Differentiating 
these equations, we get several more linear equations. Eliminating S and its derivatives from 
the system thus obtained, we arrive at an overdetermined system of non-linear ODEs for R. 
Investigating the compatibility of latter system, we find all possible functions R. Given R the 
corresponding cubic polynomial S can be easily found from (I5.66p . 

One of possible pairs R, S in Case A is given by Case A-l: 
where S3 and W are arbitrary polynomials of degree 3 and 2 (cf. Example 1). 
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It is possible to show that otherwise we have Case A-2: R = 1/W up to transformations 
of the form (15.641) . Notice that the polynomial W can be reduced by a linear transformation 
(15.641) to one of the following canonical forms: W(x) = x(x — 1), W(x) = x 2 , W(x) = x, or 
W{x) = 1. Using ( 15.661) . we find the following solutions: 

R{ x ) = , 7T > S(x) = c\ + c 2 x(x - 1), 
x[x — 1) 

R(x) = S(x) = ci + c 2 x 2 , 
x z 

R(x) = — , S(x) — c\ + c 2 x + C3X 2 . 

x 

For W = 1 see Case B. 

It is easy to verify that in Case B the only solution of (I5.65p . (l5.66p is R = 1, S = S/i, where 
5*4 is an arbitrary fourth degree polynomial. 

Now we should find the functions A and B for all above cases. It can be verified that A and 
B satisfy the same fifth order linear equation whose coefficients are differential polynomials in 
R, S. Hence all possible functions A{u) form a vector space V of dimension < 5. According to 
Remark 7 all possible functions B belong to V and dim V has to be not less then 3. 

It turns out that in Case B we have 

R=l, S = #4, A = — -, B = — i, 

04 04 

where 5*4, A4, -B4 are arbitrary polynomials of fourth degree. In this case dimU = 5. 
In Case A-2 we obtain that dim V = 4. If deg W = 2 then a basis of V is given by 

A - A 1 A 1 A U 

A t = 1, A 2 = — , A 3 = — , A 2 



W SW' z sw 

In the case deg W < 2 a basis is -^77, i = 0, 1, 2, 3. 

For Case A-l (if the pair R, S does not belong to Cases B and A-2), then the vector space 

u l 

V is three-dimensional with a basis — , i = 0, 1, 2 (cf. Example 1). 

The following example has a structure more complicated then (15. 63ft : 
Example 4. 

(u - a) 2 (z -a)(z- b) ((a + b)uz - 2b 2 z - labu + b 2 (a + b) 

X(u,z) = r . 

[u — z)u 

a basis of V is u, 1, u' 1 . 
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